Environment- induced entanglement with a single photon 
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We propose an all-optical setup, which couples different degrees of freedom of a single photon, to 
investigate entanglement generation by a common environment. The two qubits are represented by 
the photon polarization and Hermite-Gauss transverse modes, while the environment corresponds to 
the photon path. For an initially two-qubit separable state, the increase of entanglement is analyzed, 
as the probability of an environment-induced transition ranges from zero to one. An entanglement 
witness that is invariant throughout the evolution of the system yields a direct measurement of the 
concurrence of the two-qubit state. 

PACS numbers: 



I. INTRODUCTION 

The investigation of the impact of the environment on 
multi-partite entangled systems has brought out some 
subtle aspects of entanglement. Its time-dependent be- 
havior can be much different from that of the decay of 
the populations of the individual parts, or of the coher- 
ence between them. In fact, it has been demonstrated 
theoretically 0, 0, ijl MMM EL S M EH, GJ, El Q 
and experimentally [la . Il6l . Il7l ] that entanglement may 
vanish much before coherence disappears. 

Linear optical devices have been shown to be a pow- 
erful tool for the investigation of this phenomenon (Tol . 
[13, EH ■ By associating the polarization degrees of free- 
dom with qubits, and the momentum degrees of freedom 
with the environment, one has been able to realize exper- 
iments that probe the dynamics of entanglement of qubit 
systems interacting independently with individual envi- 
ronments, which, as an added bonus, can be tailored at 
will and continuously monitored. In particular, "entan- 
glement sudden death" [T^] was demonstrated [HI, and 
the evolution of a qubit under continuous monitoring of 
the environment was investigated [i~7l ]. 

The assumption of independent environments requires 
that the individual systems be sufficiently separated; for 
atoms interacting with the reservoir of electromagnetic 
field modes, this would mean that they should be much 
farther apart than the wavelength of the resonant radia- 
tion. If this condition is not met, the interaction of the 
parts with a common environment may give rise to col- 
lective effects. In this case, the direct interaction between 
the individual systems may substantially affect the collec- 
tive behavior of the system, as has been extensively dis- 
cussed within the realm of superradiance [H E(|: thus, 
for instance, when atoms get closer together, dipole in- 
teractions should be taken into account. 

In some situations, however, the direct interaction be- 
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tween the individual components of the system does not 
mask the effects of the interaction with a common en- 
vironment. This is the case, for instance, in ion trap 
experiments, where neighboring ions may be subject to 
a common dephasing process. Under these conditions, 
it is possible to protect the system against decoherence 
by defi ning log ical qubits in decoherence- free subspaces 
[2l|, [H, HI El, El, E3, El, El , which have been ex- 
perimentally demonstrated with twin-photon beams [3(| , 
nuclear magnetic resonance [3ll |32| , and trapped ions 

[11 El EE IS. 

A common environment may ent ang le initially sepa- 
rable states p3, El EE EE EH, El El H Ei EE S3, 

even when it is in a thermal state [371 . |38| . Two problems 
arise however in possible experimental demonstrations of 
this property: this effect can be masked by the direct 
interaction between the individual parts of the system 
and it may be difficult to prepare the qubits in arbitrary 
individual states due to their proximity. 

Here we show that an all-optical setup that combines 
different degrees of freedom of a single photon can be used 
to investigate the role of a common environment in the 
generation of entanglement between two non-interacting 
qubits, which can be prepared in independent individual 
states at will. Polarization and transverse mode degrees 
of freedom of the photon stand for the two qubits, and 
the photon path is the environment. Coupling between 
the qubits and the environment is achieved with linear 
optics devices. 

The use of multiple degrees of freedom of photons, 
while not leading to scalable quantum computation 
[H[, has allowed the study of basic quantum algo- 
rithms [IE E(| EH, quan tum teleportation [52|, entan- 
glement purification [53j , improved Bell-state analysis 
[HI, [55L E3] , creation of high-dimensional entangle- 
ment |5a|. demonstration of direct methods for measur- 
ing entanglement [53. | 60l|, e ngineering of mixed states 
through decoherence [la , [6l| , quantum key distribution 
[62l . l63j , measurement of a topological phase 64 1 and the 
investigation of the environment-induced dynamics of en- 
tangled systems [TH, E3 • 
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Our proposed setup is divided into two parts. The first 
one describes the action of a common environment act- 
ing on the two qubits; the second analyses the generated 
state through the measurement of an entanglement wit- 
ness, which has a peculiar and very useful property: it is 
time- independent, so the same setup can be used, as the 
environment-induced transition probability ranges from 
zero (corresponding to the initial time) to one. 

This paper is organized as follows. In Section [IT] we 
present the master equation that describes the inter- 
action of an ensemble of qubits with a common envi- 
ronment. The corresponding unitary map, essential for 
building the experiment, is introduced in Section IlIII In 
Section [TV] we present the single observable that allows 
for a direct measurement of entanglement, thus avoid- 
ing the cumbersome process of tomography [65|. The 
proposed experiment, using a linear optical setup, which 
leads to entanglement creation, plus the measurement 
circuit to quantify it, are presented in Section [Vj In Sec- 
tion lVII we show that the same optical circuit can be used 
to probe the coherent interaction between two two-level 
systems and a single-mode environment. The conclusions 
are presented in Section [VIII In the Appendix we derive 
the Kraus operators from the master equation describing 
the interaction of two qubits with a common environ- 
ment. 



II. THEORETICAL MODEL 

We consider a system of N identical qubits interact- 
ing with a common zero-temperature environment. The 
state of each qubit is represented in the basis {|e), \g}}, 
where we take |e) to be the excited state and \g) the 
ground state. The state |e) decays into the state \g), pro- 
ducing one excitation in the environment. We may asso- 
ciate to qubit i the Pauli operator Sf = \ei)(ei\ — \gi)(gi\ 
and the ladder operators = \ei){gi\ and S~ = \gi){ei\. 
The collective interaction between the qubits and the in- 
finite modes of a bosonic reservoir is described in the 
rotating-wave approximation by the Hamiltonian 

N 

H = y HuS* + ^ hue (a\a e + 1/2 



N 



--S~a\ 



(1) 



where ge is the coupling constant (taken to be real), w 
is the transition frequency between the two qubit levels, 
ai and at are the annihilation and creation operators 
corresponding to mode I with frequency uif,. Here we 
have neglected the spatial dependence of the coupling. 

The interaction of the qubits with the common envi- 
ronment is mediated by the collective operators = 

N 

Sf 1 . The dynamical evolution of the qubit system 



can be described by a master equation, which in the in- 
teraction picture, under the usual Born-Markov approxi- 
mation, and neglecting the direct interaction between the 
qubits, is given by [2dl |66| : 



Ps (t)S + S- - S+S- Ps (t)], 



(2) 



where ps (t) is the density matrix of the system and the 
spontaneous decay rate T is proportional to the square 
of the coupling constant evaluated at the transition fre- 
quency to. Thus, for instance, for a two-level atom in- 
teracting with the electromagnetic field, one may take 
I = {k, s}, where k is the wave vector corresponding to a 

plane-wave mode and s is the corresponding polarization 

i 

index. Then, gi = (uJk/2eohV) 2 eV • ft, where V is the 
quantization volume, e^ s is a polarization vector, and p is 
the transition dipole moment. Then T = oj 3 p? /SitcqHc 3 . 
In the following, we specialize our discussion to two-qubit 
systems. 



III. UNITARY MAP AND KRAUS OPERATORS 

The interaction of a system S with an environment £ 
can be described in terms of a unitary evolution encom- 
passing the system and the environment together: 



(3) 



where \<fi l )s (i — 1...4) are orthogonal states of the sys- 
tem, \p)s are orthogonal states of the environment, and 
M M are the so-called Kraus operators , acting only on 
the states of the system S, and satisfying 



(4) 



As stressed in [T7| , the map ([3]) is actually more general 
than the master equation ([2]), which is obtained from 
it, when the environment has many degrees of freedom, 
under Markovian and differentiability assumptions (68l . 

The Kraus operators corresponding to Eq. @ are ob- 
tained in the Appendix, using the Choi matrix formal- 
ism [2O, [7l|. Substituting the resulting Kraus opera- 
tors (|35H38p in Eq. ([3]), we get, in the collective basis 
{|0,0) ee _i=(|e 5 )-| ffe )),|l,l) ee |ee),|l,0) ee -^(\eg) + 

\ ge )),\l,-l) = \gg)}: 

|l,l)s|0) fi A|l,l) s |0)e + fl|l,0)s|U>£ + 
D\l,0) s \l B ) £ +F\l,-l) s \2) £ , 

|l,0)s|0) f - i4|l,0)s|0> £ + C|1 I -1) S |1 A )£ + 
E\l,-l) s \l B ) £ , 
|l,-l>s|0> £ - |1,-1> S |0) £> 

|0,0) s |0) £ |0,0) 5 |0) £ , (5) 
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where A,B,---F are time-dependent coefficients explic- 
itly calculated in the Appendix and the states of the 
environment arc labelled so as to highlight their phys- 
ical meaning in terms of the number of excitations, i.e., 
\0}s corresponds to the vacuum state, \Ia)s and |ls)g to 
states with just a single excitation and \2)s to a state with 
two excitations in the environment. This set of equations 
has an immediate physical interpretation: each coeffi- 
cient represents the probability amplitude of transition 
between the qubit states as a function of time, either the 
qubits emitting excitations into the environment, which 
is initially in the vacuum state, or exchanging excitations 
between them, or yet not emitting any excitation at all, 
when both system and environment remain in the same 
state. For instance, in the first line, where both qubits 
are initially in the excited state |1, 1)5, we have a proba- 
bility amplitude A that both qubits remain in the excited 
state, with no excitation being emitted into the environ- 
ment, which remains in the vacuum state |0)g. On the 
other hand, there is a probability amplitude B(D) that 
a single qubit decays, the state of the system going to a 
superposition of \eg)s and \ge)s with just one excitation 
being transferred to the environment, which goes to state 
\^a)s Finally there is a probability amplitude F 

that both qubits decay to |1, — 1)5 emitting two photons 
into the environment, whose state becomes \2) £ . Similar 
reasonings can be applied to the other lines of the map. 

We note, in Eq. ([5]), the appearance of two orthogonal 
environment states with just a single excitation, \\a) £ 
and |1_b)£. This can be understood physically by rewrit- 
ing Eq. ([5]) as: 

|M)s|0)£ -> A|l,l> s |0) £ + G|l > 0> s |li,i>£ + 
F|l,-l) s |2) s , 

|l,0) s |0) fi A|l,0} s |0}£ + if|l,-l>s|li,o>£, 
|l,-l>s|0>£ -> |1,-1) 5 |0) £ , 
|0,0) s |0) £ - |0,0) s |0) £ , (6) 

where |li i)g and 1 1 1.0) e arc non-orthogonal states of the 
environment with a single excitation given by: 



II 



B D 

q\1a)£ + ^|ls>£ 



G = VB 2 + D 2 

|ll,0>£ = 



C E 
Jj\1a)s + Jj\Ib) £ ) , 



Calculating this rate for the initial states |1,1) and 
1 1, 0) we get: 



wl' 1 = r(i,i|5+5 _ |i,i) =2r, 

Wl'° = Y{l,Q\S + S-\l,Q) = 2T. 



(8) 



The rate corresponding to |1, 1) is exactly the one we 
would expect from two independent qubits. This implies 
that in this case the first photon emitted would have 
a lincwitdh of T. The decay rate of the single excited 
qubit in the state |1,0) is twice that of each qubit in 
|1,1), so this state is superradiant and the linewidth of 
the emitted photon is expected to be 2r. As we have 
two different linewidths for the photons emitted from the 
initial states |1, 1) and |1,0), this explains why we have 
two different environment states corresponding to a single 
excitation, |li i)g and |li,o)f m the first and the second 
line of Eq. ©, respectively. 

In terms of the computational basis 
{|ee), I eg), \ge), \gg)}, Eq. ([5]) can be rewritten as 



|ee) <s |0) £ - 

\eg)s\0)s - 
l<?e) 5 |0) £ - 

\gg}s\o} £ - 

where M 



M\ee) s \0) £ + P{\eg) + |<?e)) s |l ee ) f + 
N\gg) s \2) s , 

Q\eg) s \Q) £ + R\ge) s \Q) £ + S\gg) s \l eg ) £ 
Q\ge) s \Q) £ + R\eg) s \Q) £ + S\gg) s \l eg ) e 
\99)s\0)s, 



(9) 



-rt 



P 



-2Tt 



2Yte 



-2Tt 



Vrte- 2rt , 



N = 



l_ e -2rt 



. The new single-excitation states |l ee )f an d 
\l e g) £ are related to orthogonal states of the environment 
|1a>£ and \l B )e by: 



Lee/£ 



P 

5' 



LA>£ 



AS 



w 
5" 



1 



B) £ 



(10) 



H = VC 2 + E 2 . 



where X, Y, Z and W are time-dependent coefficients 
defined in the Appendix. 

As shown in Section [V] Eq. ([9]) leads to a simple linear 
optics demonstration of environment-induced entangle- 
ment. 



IV. ENTANGLEMENT GENERATION AND 
MEASUREMENT 



The explanation for the presence of two different states 
with a single excitation in the first and second lines of 
Eq. ([6|) stems from the consideration of the photon emis- 
sion process for the initial states |1, 1) and |1, 0), respec- 
tively. The rate of photon emission for the two qubits is 
given by [20j | : 

w 2 = r(s+s-). (7) 



A detailed analysis of the generation of entangle- 
ment between two atoms via spontaneous emission, with 
dipole-dipole interaction, was presented by Ficek and 
Tanas [40 . [47j . Here we consider the special case in which 
the separation between the atoms is much smaller than 
the typical wavelength of the emitted radiation, and there 
is no dipole-dipole interaction. Our expressions coincide 
with those in Refs. 0, [4^ in the proper limit. 
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For the quantification of entanglement we use the con- 
currence [721 ] . defined as: 

C = max{0,A}, (ff) 

where: 

A= VAT- v 7 ^- Vh- y/Xi, (12) 

Xi being the eigenvalues in decreasing order of the matrix: 

p(<r y ®<r v )p*(<T y ®cry), (13) 

where p is the two-qubit density matrix, a y is the second 
Pauli matrix and the conjugation is performed in the 
computational basis. Concurrence ranges from 0, which 
corresponds to a separable state, to 1 which corresponds 
to a maximally entangled state. 

If the system is initially in the separable state ps(0) = 
\eg)(eg\, then, by applying Eq. |9]) and tracing over the 
states of the environment, we end up with the following 
state for the system: 

/e- rt + l\ 2 /p- rt -l\ 2 

ps{t)= [-^^) \eg){eg\ + [— r - J \ge)(ge\ 

e~ 2Tt - 1 
+ 1 (\eg){ge\ + \ge)(eg\) 

+ 1 ~ e 2 2Vt \gg)(gg\- (u) 

Simple calculation shows that the concurrence for this 
state is given by 

C(t) = i(l- e - 2rt ). (15) 

This clearly shows that, although initially in a separa- 
ble state, the system evolves to an entangled state, its 
concurrence reaching the maximum value of 1/2 in the 
asymptotic regime t — > oo. This entanglement is induced 
solely by an indirect interaction mediated by the common 
environment, since in the model considered here there is 
no direct interaction between the two qubits. Therefore, 
in order to observe the creation of entanglement it suf- 
fices to implement experimentally the second (or third) 
line of Eq. ©. 

The emergence of entanglement in this case can be 
qualitatively understood by considering that the initial 
separable state \eg) can be expressed as a superposition 
of the singlet |0,0) and triplet |1,0) Bell states [73| : 

|e 5 ) = -^(|0,0) + |l,0)). (16) 

One can trivially see from Eq. ([2j that the singlet state 
does not evolve. On the other hand, the triplet state 
decays asymptotically to the ground state. So as time 
approaches to infinity we end up with the following mixed 
state: 

P«t = i|l,-l><l,-l| + i|0,0><0,0|. (17) 



The concurrence of the above state is easily calculated 
to be 1/2, in accordance with the asymptotic limit of 
Eq. (|15j) , the contribution for entanglement coming from 
the maximally-entangled singlet state in Eq. {T7|. 

One should note that for two closely-spaced atoms, the 
initial state \eg) is usually difficult to prepare, since it 
does not have the same symmetry as the ground state 
under exchange of particles, and the interactions with 
external fields are symmetrical, as the interaction with 
the environment. 

On the other hand, if both qubits are initially excited, 
i.e, ps(0) = |ee)(ee|, which is in principle simpler to 
prepare, one should not expect to have asymptotic en- 
tanglement, since this initial state does not have a sin- 
glet component. The question remains however whether 
transient entanglement might still occur. At intermedi- 
ate times, according to the first line of Eq. ([6|), the state 
of the system becomes a mixture of the initial state, the 
maximally- entangled state (\eg) + \ge))/y/2, and the state 
\gg). In spite of the presence of the maximally- entangled 
component (which eventually decays to | <?<?)), the state 
remains separable for all times. Indeed, the concurrence 
of ps{t) is given by 

C(t) =max{0,Ci,C 2 }, (18) 

with 

d{t) = -2e- 2rt Tt, 

Ca(t) = 2e- 2rt (Tt- \/e 2rt -2Tt-l). (19) 

As Ci(t),C2(t) < 0, no entanglement is generated in this 
case. The same conclusion is obtained from Ref. [4(| 
when the distance between the atoms goes to zero - 
as shown in that reference, the dipole-dipole interaction 
does not play any role for this kind of initial state. 

We will thus concentrate on the state described by 
Eq. ([14]) . Even though this state is difficult to prepare 
in most systems, as commented above, this is not the 
case in the experiment proposed here, as will be shown 
in the next Section. Using this state has an important 
and useful consequence. Instead of doing tomography 
[6l| for full reconstruction of the state, in order to deter- 
mine the entanglement through Eq. (|11[) . which would 
require measuring sixteen observables, we can define a 
time-independent entanglement witness - a unique ob- 
servable that could be used to detect entanglement pro- 
vided that we have some a priori knowledge about the 
states involved. 

By definition, an entaglement witness W is an observ- 
able that satisfies the following relations: Tr(W p sep ) > 
for any separable state and Tr(Wp en t) < for at least 
one entangled state, which is then said to be detected by 
the witness. 

For an arbitrary bipartite density matrix p. an entan- 
glement witness is given by (74[: 

W = 1 - {UV^) Ta , (20) 
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where U and V are unitary matrices obtained from the 
singular value decomposition of p TA , Ta being the partial 
transposition with respect to any one of the two qubits. 
The singular decomposition is given by p TA = UYiV^ , 
where £ is a diagonal matrix whose elements are the 
nonncgative square roots of the eigenvalues of p TA p T - A ^ 

m- 

For the family of states given by Eq. (fl~4|) , evaluation of 
Eq. (|2"0|) leads in fact to a time-independent entanglement 
witness: 



W = 



-^\9e){eg\ 



1 



71 ' 



(21) 



Furthermore, this witness is optimal for the familiy 
of states here considered since it leads directly to the 
concurrence, which can be calculated straightforwardly 
in terms of the mean value of W as: 



C = 



Tv(W Ps (t)) 
(1-V2) 



(22) 



Therefore, by measuring just a single observable we are 
able in this cast 1 not only to detect but also to quantify 
entanglement, throughout the evolution of the system. 



V. EXPERIMENTAL SETUP 

We discuss now the experimental setup, involving 
linear optics, which leads to the demonstration of 
environment-induced creation of entanglement, and also 
the setup for measuring the entanglement witness in Eq. 
(|2Tj) . It is possible to design a circuit that implements the 
whole set of equations in Eq. ([9]), but for our purposes it 
is enough to consider just the second line of Eq. (|9|). 

Our strategy for studying the dynamics of entangle- 
ment is similar to the one adopted in Refs. (TBI . Il7j |. 
There, the unitary quantum map corresponding to the 
evolution of a single decaying atom was described in 
terms of the decay probability p, which is an exponential 
function of time. This allowed us to implement the de- 
cay dynamics through a static linear optics setup with a 
changeable parameter, namely the orientation 9 of a half- 
wave plate, with p = sin 2 (26). It also allows, from the 
theoretical point-of-view, to encompass different kinds of 
dynamics in the same map, since p could also be consid- 
ered to be an oscillating function of time, and in this case 
the quantum map would correspond to Rabi oscillations 
of a two-level system interacting with a single mode. 

Here we generalize this approach by parametrizing the 
evolution in terms of two fundamental probabilities: the 
probability that the first qubit decays emitting one ex- 
citation into the environment; and the probability that 
the excitation goes from the environment to the second 
qubit. Different values of the coupling constant in Eq. ([TJ) 



TABLE I: System encoding using polarization and Hermite- 
Gaussian modes. 



|ee) = 


\V,HG 10 ) = 


\Vv) 


\eg) = 


\V,HG i) = 


\Vh) 


\ge) = 


H, HG\o) = 


\Hv) 


\99) = 1 


H, HGm) = 


\Hh) 



result, for the same instant of time, in different values of 
these probabilities, which can be tuned by changing an- 
gles 6*i of a half-wave plate and 62 of a Dove prism in the 
proposed experiment. Alternatively, one may say that 
different values of Ft correspond to different values of 0\ 
and 62, since the coupling constant always shows up in 
this combination. Of course, these two angles are not 
independent, since they are related through their respec- 
tive dependences on Tt. 



A. SYSTEM-ENVIRONMENT EVOLUTION 
CIRCUIT 

Both the system and the environment are represented 
by different internal degrees of freedom of a single photon. 
We associate the horizontal (H) and vertical (V) polar- 
ization of the photon to the ground and excited state of 
the first qubit respectively and the first-order Hcrmite- 
Gaussian transverse modes [76[ to the states of the second 
qubit - mode HGqi is associated to the ground state and 
HG 10 to the excited state. For simplicity, we define the 
modes as: HG01 = \h) and HGm = \v). The environ- 
ment states arc represented by different paths of the pho- 
ton. The encoding for the system is summarized in Table 
HJ Figure Q] shows the proposed experimental setup. It 
employs optical components that allow for independent 
manipulation of polarization, transverse and longitudi- 
nal (path) spatial degrees of freedom. Since in this case 
one has only one environment state with one excitation, 
\leg)s is redefined as \l)s- The circuit is composed of a 
spontaneous parametric down-conversion source, a holo- 
graphic mask (HG Mask), a half- wave plate (HWP1), 
a Mach-Zender interferometer with an additional mirror 
(MZIM1) [z3, a Dove prism (DPI), two polarizing beam 
splitters (PBS1 and PBS2) and two glass plates (GP1 
and GP2). 

The MZIM is shown in Fig. [2j For an optical phase 
difference between the arms properly adjusted by tilting 
a glass plate (GP), \Vh) and \Hv) (\Vv) and \Hh)) states 
entering the MZIM through port 1 leave it from port 4(3) 
and the \Vh) and \Hv) states (\Vv) and \Hh)) entering 
through port 2 exit from port 4(3). 

The parametric down-conversion source, composed by 
a non-linear type-I cristal (Cr) pumped by a horizontally- 
polarized laser beam, produces two vertically-polarized 
twin photons. The two-qubit system plus environment 
is encoded in the lower photon. The upper photon will 
solely be used to validate photon counts in the lower path 
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FIG. 1: (Color online) Optical setup to investigate environment- 
induced entanglement creation. HG Mask stands for holographic 
mask , HWP for half-wave plate, MZIM for Mach-Zender inter- 
ferometer with an additional mirror, DP for Dove prism, PBS for 
polarizing beam splitter, GP for glass plate, M for mirror, and D 
for photon detector. 



The physical interpretation of the transformation l|24[). 
performed by the dashed circuit block A, is that there 
is a probability amplitude cos(2#i) that the system and 
environment remain in the same state | e^) s 1 0) £ , and a 
probability amplitude sin(2#i) that the first qubit decays 
emitting one excitation into the environment. 

In the path associated to the state of the environment 
with one excitation \l)s, the photon propagates through 
a Dove prism rotated at an angle 9 2 , that accomplishes 
the following transformation: 



Infill) 



[cos(28 2 )\Hh) + sm(20 2 )\Hv)] s \l)e- (25) 



After the Dove prism DPI, the Mach-Zender interferom- 
eter MZIM1 conserves the path of the \Hh) state and 
changes by 90° the path of the \Hv) state. We asso- 
ciate the horizontal path to the vacuum state of the en- 
vironment which we call \Q")s and the vertical path after 
the MZIM1 to the environment state with one excitation 
The resulting transformation after the MZIM1 is 
the following: 



\Hh) s \l)s -» cos(20 2 )|fl7i) 5 |l) £ + sin{26 2 )\Hv) s \Q") £ . 



(26) 




FIG. 2: (Color online) Mach-Zehnder interferometer with an ad- 
ditional mirror. BS stands for beam splitter, GP for glass plate, 
and M for mirror. 



in a coincidence detection regime. The glass plate GP1 
is used to adjust the upper and lower optical paths. To 
generate the initial state \eg) = \Vh) we place the holo- 
graphic mask, specially designed to produce the HGqi 
mode, in the path of the lower photon. The propaga- 
tion direction of the photon after the mask is associated 
to the vacuum state of the environment |0)f. After ini- 
tial preparation in state |V7i)s|0)£ the photon propagates 
through the half- wave plate HWP1 aligned at an angle 
6\ with respect to the vertical polarization. The trans- 
formation perfomed by HWP1 is: 

\Vh) s \0)e -> [cos(2fli)|Vfc) + sin(20i)|JTfc}] s |0) s . (23) 

After HWP1, the polarizing beam splitter PBS1 reflects 
(transmits) photons in the \Vh) (\Hh)) state. We asso- 
ciate the horizontal path after PBS1 to the vacuum state 
of the environment, which we call |0')g, and the vertical 
path to the environment state with one excitation 
This leads to: 



\Vh) s \0)e -> cos(20i)|Wi> s |O')£ + sin(20i)|IT&) s |l) £ . 



(24) 



The physical interpretation of the transformation (|26| . 
performed by the dashed circuit block B, is that there 
is a probability amplitude cos(26>2) that the system and 
environment remain in the same state, corresponding to 
the two qubits in the state \g) and one excitation in the 
environment, and a probability amplitude sin(2f52) that 
the excitation goes from the environment to the second 
qubit. 

Finally, the polarizing beam splitter PBS2 together 
with the GP2 allows for the in-phase coherent combina- 
tion of the two paths associated with the vacuum state 
of the environment |0')f and \0")e into a single horizon- 
tal path |0)g as shown in Figure [T] Therefore the overall 
transformation of the whole circuit is: 

\Vh) s \0) £ ^ [coa{2e 1 )\Vh) 

+ sin(20 1 )sin(20 2 )|i^)] 5 |O) £ 

+ sin(20i) cos(26 2 )\Hh) s \l)s ■ (27) 

This equation coincides precisely with the second line 
of Eq. ©, upon associating cos(26*i), sin(26*i) sin(26>2) 
and sin(2#i) cos(2f?2) with the coefficients Q, R, and S in 
that equation, respectively. 

Note that the coefficient Q is the probability ampli- 
tude that the first qubit remains in the excited state, R 
is the probability amplitude that the first qubit decays, 
emitting one excitation into the environment, multiplied 
by the probability amplitude that the second qubit ab- 
sorbs the excitation present in the environment, and S 
is the probability amplitude that the first qubit decays 
multiplied by the probability amplitude that the second 
qubit remains unexcited, so that the excitation emitted 
by the first qubit remains in the environment. 
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The interaction of the system with the environment 
is thus achieved in a entirely controllable manner: each 
value of Tt corresponds to a specific set of angles 9\ and 
6*2- Furthermore, this setup makes it clear that the inter- 
action between the two qubits is mediated by the envi- 
ronment. One should also note that the present proposal 
allows for the easy preparation of the initial state \eg), 
which, as mentioned before, may be a demanding task in 
other systems. 



path). Rotated at 45° the Dove prism in the upper 
path flips between Hermite-Gaussian modes: HGqi <-> 
HGio = h <-> v. The optical lengths of the upper and 
lower paths are adjusted by the glass plate GP3 enabling 
coherent in-phase combination in PBS4. This arrange- 
ment corresponds to the usual CNOT gate where the po- 
larization degree of freedom acts as the control bit and 
the transverse mode acts as the target bit, according to 
Table HI 



B. MEASUREMENT CIRCUIT 

In order to measure the mean value of the entangle- 
ment witness (|2"Tj) we decompose it in its diagonal basis, 
which happens to be the collective basis: 



W = ( 1 + 7l) IMXMI + ^IMXMI- 
-L|o,o)(o,o| + (i-^=) |i,-i><i,-i] 



TABLE II: CNOT gate. C stands for control and T for target. 



(28) 



In Fig. [3] we show the circuit that performs projective 
measurements of ps{t) in the collective basis. Photons 
in modes \0}s and outgoing the system-environment 
evolution circuit are incoherently driven to the measure- 
ment circuit, which corresponds to a trace over the states 
of the environment. 

The circuit is composed of a Mach-Zender interferom- 
eter with an additional mirror (MZIM2), two polarizing 
beam splitters (PBS5 and PBS6), a CNOT gate 0, a 
half-wave plate (HWP2), and four photon detectors (Dl, 
D2, D3 and D4). 



MZIM2 



|0>._ 
|1>." 



PBS6 



PBS3___pP2@45°GP3 PBS4 H WP2@22.5 PBS5 D1 



CNOT 



D2 



D3 



I D4 



FIG. 3: (Color online) Optical setup to measure the entanglement 
witness. MZIM stands for Mach-Zender interferometer with an 
additional mirror, PBS for polarizing beam splitter, DP for Dove 
prism, HWP for half-wave plate, GP for glass plate, M for mirror, 
and D for photon detector. 

The CNOT gate is an interferometer composed by 
two polarizing beam splitters (PBS3 and PBS4), a Dove 
Prism (DP2) and a glass plate (GP3). PBS3 splits hori- 
zontal (upper path) and vertical polarized photons (lower 



Input 
C T 
\Vv) 
\Vh) 
\Hv) 
\Hh) 



Output 
C T 
\Vv) 
\Vh) 
\Hh) 
\Hv) 



The MZIM2 conserves the path of the \Vh) and \Hv) 
states and changes by 90° the path of the \Vv) and \Hh) 
states. This implies that photons in states |1, 0) and |0, 0) 
will be transmitted by the MZIM2 while photons in states 
|1, 1) and |1, —1) will be reflected by it. 

In the horizontal path after the MZIM2, the CNOT 
gate circumscribed in the dashed block and described 
by Table HH plus HWP 2 aligned at 22.5° with respect 
to the horizontal polarization, which plays the role of 
a Hadamard gate in the polarization degree of freedom, 
i.e., \H) -> -fe(\H) + \V)) and \V) - j-{\H) - \V), per- 
forms the following transformations: |1,0) — > \Hh) and 
|0,0) -> -\Vh). Then, by virtue of PBS5, which trans- 
mits the horizontal polarization and reflects the vertical 
one, Dl and D2 detect the populations of |1, 0) and |0, 0), 
respectively. In the vertical path after the MZIM2, due 
to PBS6, populations of |1, 1) and |1, —1) are detected in 
D3 and D4, respectively. 

Let Ci be the photon counts in detector i and N = 
2j=i Ci the total counts during the measurement inter- 
val. Then, according to Eq. 



Tr(W Ps (t)) 



C 2 

+C 4 ( 



C 3 ll 



i 



1 



(29) 



The concurrence follows straightforwardly through Eq. 



VI. SINGLE-MODE ENVIRONMENT AND 
INITIAL ENTANGLED STATES 

Our proposal can be easily generalized to describe the 
evolution of other initial states as well as the oscillatory 
exchange of energy between the two-qubit system and a 
single mode of the environment. As mentioned in Sec- 
tion [Vj the approach of unitary quantum maps is quite 
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general and accommodates other types of dynamics, de- 
pending on the way we parametrize the coefficients with 
time. Indeed, Eq. j9]) can match the non-dissipative dy- 
namics of two non-interacting qubits in a resonant cavity, 
described by the following Hamiltonian: 

2 2 
H = y ^huSf + hw((Ja+l/2)-ihgYj ({Sfa- Srat) . 

i=l i=l 

(30) 

where g is the coupling constant, while a and are the 
annihilation and creation operators of the cavity field 
mode. In this situation |l ee )£ and |l e g)£ are identified, 
since there is just a single mode in the cavity. No ma- 
jor changes, except for the time dependence of the co- 
efficients, are needed to represent this rather different 
dynamics. As a matter of fact, the circuit describing 
the dissipative dynamics can be applied to the resonant 
case as long as the coefficients of the map in Eq. (|9]) 
are given as functions of time by: M = h(l + cos (2gt)), 
P = -±wx(2gt), N = i(-l + cos(2$t)), Q = ±(1 + 
cos (y/2gt)), R = §(-1 + cos (y/2gt)), S = ^|sin {yfegt). 
The extra phase i can be easily compensated by a phase 
plate. 

The angles 9\ and 82 in the proposed experimental 
setup are now functions of gt. Tuning the coupling con- 
stant, or rather the product gt, corresponds to tuning 
those two angles, according to the above equations. 

Our setup encompasses therefore both the decay and 
the oscillatory case, as in the experiment reported in 
Ref. [l5j . This is a great advantage of our parametriza- 
tion of the evolution of the system. 

Similar circuits can be applied to analyze the dynamics 
of initially entangled states, which are easily prepared 
within the present framework. For instance, the state 
a\eg) + (3\ge) is obtained from the initial state \eg) by 
applying a half-wave plate to the polarization mode and 
the CNOT circuit in Fig. [3] in order to entangle the two 
degrees of freedom. 

The dynamics of the family of states a\gg) + (3\ee) un- 
der the action of a global environment is particularly rich, 
giving rise not only to the disappearance of entanglemcmt 
at finite times but also to entanglement revival. A de- 
tailed analysis was presented in Ref. [47| . where it was 
shown that in this case the evolution of the density ma- 
trix elements is independent of the dipole-dipole interac- 
tion between the atoms. Therefore, the results obtained 
from our model precisely coincide with those in Ref. p7| , 
in the limit when the distance between the two atoms 
goes to zero. 



VII. CONCLUSIONS 

In this work we propose a simple experiment, which 
uses only elementary linear optical devices, to observe 
creation of entanglement between two non-interacting 
qubits, induced by a common environment. 



We encode the two qubits and the environment into 
different degrees of freedom of a single photon: the po- 
larization, the transverse modes, and the different paths 
stemming from different directions of the photon momen- 
tum. Our experimental setup directly implements the 
unitary map that describes the time evolution of the sys- 
tem and the environment. We recover the dynamics of 
the two-qubit system by tracing over the path of the pho- 
ton, which represents the environment. 

The evolution map is parametrized in terms of time- 
dependent transition probabilities, rather than time it- 
self. This allows the investigation of the dynamics of 
entanglement through a static setup. It also follows from 
this parametrization that the same setup describes not 
only the collective interaction of two qubits with a zero- 
temperature bath, but also the reversible exchange of ex- 
citation between the two qubits and a single cavity mode. 

Entanglement is quantified by a time-independent en- 
tanglement witness - a single observable - shown to 
be proportional to the concurrence for the states here 
considered. Therefore, our proposed scheme to quan- 
tify entanglement does not require time-consuming to- 
mographic measurements to determine the state of the 
system. 



VIII. APPENDIX 

A. Kraus Operators 

The first step to obtain the Kraus operators is to deter- 
mine ps (t) solving the Lindblad equation @ . In the col- 
lective basis {|0, 0) = -±(\eg)-\ge)),\l,l) = |ee>, |1,0> = 
7g(l e 0) + l^e)), 1 1, — 1> = \gg)}, which we denote by the 
indices 1,2,3,4, respectively, ps(t) has the following ele- 
ments: 

Pn(t) = MO), 

P12® = Pi 2 (0)e- rt , 

Pis(t) = Pi 3 (0)e- rt , 

Pii(t) = pii(0), 

P22(t) = P22(0)e' 2Ft , 

P2 3 (t) = P23(0)e~ 2r \ 

P2i{t) = p 24 (0)e- r *, 

P33(t) = P3 3 (0)e- 2r * 

P34(t) = p 34 (0)e- rt - 

p i4 (t) = P44(0)+/922(0)(l-e 

+P3 3 (0)(l-e- 2r *). (31) 



-2rt 



2p 22 (0)Tte 
2p 23 (0) e - r *(l-e- r *) 



-2rt 



2rte 



-2rt\ 



We define now the matrix E l i (t) as ps(t) with the 
initial conditions p m „(0) = S m iS n j. For example, E 22 (t) 
is given by ps(t) considering solely /?22(0) = 1 and all 
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other coefficients equal to zero. Using (f3Tj) we get: 



E 2 \t) 



/o o o\ 

e" 2rt 

2Tte- 2rt 

\ ) 



(32) 



The next step is to build the 16 x 16 positive scmidcf- 
inite Choi matrix [7l[ composed of matrices E l] (i) as 
follows: 



E u (t) E 12 (t) 
c= I E*\t) E**(t) 



Matrix C can be represented by: 



(33) 



where: 

A - ~- rt 



e -,B = C 



VI71 ; 



H2 JJ _ Pi P2 _ ft 

x/2 7l ' v^72 ' %/272 



F = V I - e- 2rt - 2e- 2rt - Tt, 

O = ^17 - 32e rt + e 4r * + e 2r *(14 - 4Tt) + 4Ti(l + R), 



7i 

72 



i-e Jit +2rt-n 

4(-l+e rt ) ' 

l-e 2rt +2rt+Q 
4(-l+e rt ) ' 



= v /e- 2rt (-l + e 2rt + 2R-fi), 
= v / e -2rt(_ 1 + e 2rt + 2 n + fi), 

-, . (~l+e 2 "-2rt+Sl) 2 
1 16(-l+e rt ) 2 ' 



(l-e 2rf +2rt+n) 2 



16(-l+e rt ) 2 



(34) 



where 



are the normalized eigenvectors of C with 



eigenvalues a^, and are non- normalized vectors. In 
the present case we have just four non-null vectors {v^}, 
each one corresponding to a Kraus operator. For conve- 
nience we designate [i = 0, 1a, 1b, 2. 

The final step is to divide each column matrix corre- 
sponding to \Vfj,) into four segments of equal length. Then 
the matrix representing Kraus operator M M has the ith 
segment of as its ith column. 

Following this procedure, we end up with the following 
Kraus operators: 



M 



( 1 

A 

A 

\ 1 



(35) 



B. Numeric Coefficients 



Here we show the time-dependent coefficients appear- 
ing in Eq. (fTU]) : 



ia. 

2£ 2 : 



-1+e 1 '*) 



Y 



X 

w 

5i = ^e- 2rt (-l + e 2rt + 2Yt - Q), 



(l-e 21t +2rt+n)<5 2 y _ gi 
8(-l + e rt )e 2 ' ~~ 2ei 



= */l 



(-l+e 2rt -2rt+Q) 2 



16(-l+e r *) 2 



5 2 = ^ e -2rt(_ 1 + e 2r t + 2 n + r2)", 

„ _ /-, | (l-e 2 »+2rt+Q) 2 
2 — V ^(-l+e 1 '*) 2 ■ 



Mia 



/ 



B 

\0 C 



(36) 
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